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Using time-dependent Ginzburg-Landau theory we demonstrate that the Aharonov-Bohm (AB) 
effect, resulting from a Berry phase shift of the (macroscopic) wavefunction, is revealed through the 
dynamics of topological phase defects present in that same wavefunction. We study vortices and 
ant ivort ices on the surface of a hollow superconducting cylinder, moving on circular orbits as they 
are subjected to the force from the current flowing parallel to the cylinder axis. Due to the AB 
effect the orbit deflections, caused by a magnetic field component along the cylinder axis, become 
periodic as a function of field, leading to strong and robust resistance oscillations. 

PACS numbers: 74.78.-w, 74.25. F-, 74.20.De 



I. INTRODUCTION 



II. MODEL 



Recent advances in nanotechnology have triggered in- 
terest in experimental and theoretical studies of meso- 
scopic and nanoscale superconductivity on curvilinear 
surfaces^"—. In superconducting spherical nanoshells the 
surface curvature leads to a Magnus- Lor entz force, which 
pushes the vortices and antivortices towards the opposite 
poles of the shell. This can be considered as an effective 
pinning of vortices and antivortices at the poles, which 
strongly affects both the equilibrium distributions of vor- 
tices and their dynamics^^. The effects of surface cur- 
vature on vortex dynamics have been recently analyzed 
also for curved stripes and hollow cylinders^. An impor- 
tant aspect of hollow cylinders, which have made them 
a popular subject of study, is their doubly connected 
topology. The Aharonov-Bohm (AB) oscillations^, orig- 
inating from a shift of the geometric (Berry) phaseii of 
the charge-carrier wave function by an enclosed magnetic 
flux, were first experimentally observed in superconduct- 
ing hollow cylinder si 2 - and ringsii In the present paper, 
we analyze the vortex dynamics in superconducting hol- 
low cylinders, subjected to a magnetic field tilted with 
respect to the cylinder axis. In particular, we investigate 
the effect of the oscillating persistent current, related to 
the AB effect, on the vortex motion and the correspond- 
ing resistive state. The obtained results suggest that, in 
parallel to the critical temperature oscillation a 13 i 14 , also 
oscillations of the resistance, caused by vortex motion, 
can be used as a tool to probe the AB persistent cur- 
rents in superconducting hollow cylinders. 



The sketch of the structure under consideration is 
shown in Fig.OJa). A thin superconducting hollow cylin- 
der with thickness length L and radius R is subjected 
to an external homogeneous magnetic field Bo, which in 
general has nonzero components both along the cylinder 
axis (Bq\\) and in the perpendicular direction (Bo±). In 
the chosen cylindrical co-ordinate frame, the z-axis coin- 
cides with the cylinder axis, while the direction of Bq±_ 
corresponds to the angular coordinate <p = tt/2. At z = 
and z = L, where normal- met al/superconducor bound- 
ary conditions are applied, an external current with den- 
sity j e is along the z-axis. 

The vortex dynamics is described within the time- 
dependent Ginzburg-Landau (TDGL) approach. For nu- 
merical simulations we use the implementation of this 
approach described in detail in Refs. [H and [HI and ap- 
plicable when the thickness of the superconductor is 
smaller than the Ginzburg-Landau coherence length £, 
with obvious adaptations to the case of a cylindrical 
layer. Like in Ref. [HI, the relevant quantities are made 
dimensionless by expressing lengths in units of \/2£, time 
in units of 7r^/[4fc#(T c — T)] « 11. §tgl, magnetic field in 
units of <!>o/(47r£ 2 ) = H c2 /2, current density in units of 
<I>o/(2V / 27r/ioA 2 = 3VS/(2y/2)j c , and scalar potential 
in units of 2/c#(T c — T)/(7re). Here, <I>o = 7rh/e is the 
magnetic flux quantum, /i is the vacuum permeability, 
A is the penetration depth, tql is the Ginzburg-Landau 
time, H C 2 is the second critical field, and j c is the crit- 
ical (depairing) current density of a thin wire or film^. 
The results, described below, are obtained for a fixed 
value of the Ginzburg-Landau parameter (k = 0.77). 
Note that this particular choice of n is not really re- 



2 



(a) 



(b) 





I 



O-F-i — i — i— 1 I i i i 

71 71 

FIG. 1. (Color online) (a) Structure under consideration. The 
red spots at <j> — tt/2 (light blue spots at (j) — 3tt/2) schemati- 
cally show equilibrium positions of vortices (antivortices), in- 
duced by the field B ±, in the case of j e = and B \\ = 0.(b) 
Equilibrium distribution of the square modulus of the order 
parameter in the hollow cylinder with L = 75, R = 3, and 
d — » at B ± = 0.35, B \\ = 0, and j e = . (c) Same as 
in panel (b) but for j e = 0.22. The voltage drop V is calcu- 
lated between the points labeled as "l"and "2" and located 
at z = L/4 and z — 3L/4, respectively. 



strictive: with the used units, the parameter n explic- 
itly enters the equations only through the expression 
Ai(r) = (27Tft 2 ) -1 / dVj(r')|r - r'| _1 , which describes 
the vector potential induced by the currents j flowing in 
the superconductor (see Ref. [l5| for more details). For 
thin superconducting cylindrical shells under considera- 
tion {d < 1, d <C R, L), the vector potential Ai is pro- 
portional to the ratio d/n 2 , so that the results obtained 
for k = 0.77 and a given value d can be applied also to 
cylinders with larger or smaller values of k, provided that 
the quantity d/ n 2 remains the same. 



III. RESULTS AND DISCUSSION 

In the absence of an applied current, the vortices and 
antivortices^, induced in a cylindrical shell by a mag- 
netic field Bo± = 0.35, are "geometrically pinned" by a 
Magnus- Lor entz forced to <fi = tt/2 and <\> = 3tt/2, respec- 
tively [see Fig.QJb)]. When applying a relatively low cur- 
rent density j e along the cylinder in the z-direction, vor- 
tices and antivortices experience a force towards </> = tt 
where they can annihilate. Some of the vortices and an- 
tivortices do annihilate with each other, but most of them 
remain "pinned", although their equilibrium positions 
are somewhat shifted towards (j) = tt due to the Lorentz 
force proportional to the applied current [Fig.QJc)]. At a 
certain current density, the system enters the dissipative 
regime: the vortex- ant ivort ex pairs continuously recom- 
bine at (j) = tt while new pairs are nucleated at <j) = 0, 
where the applied current, summed up with the Meisser 



current, significantly suppresses the order parameter [see 
Fig. 1(c)]. 

In Fig. [51[a) the critical current density ji, correspond- 
ing to the onset of the resistive state in the middle part of 
the cylinder (between z = L/3 and z = 3L/4), is plotted 
as a function of the radius R of the cylindrical shell at 
a fixed applied magnetic field Bq±_ = 0.9, B \\ = 0. As 
seen from Fig. [21(a), the critical current density ji(R) for 
cylindrical shells is appreciably higher than that for the 
corresponding flat stripes with the same cross-section and 
manifests a nonmonotonous behavior. As the radius is 
decreased (down to R « 1.5), the critical current density 
increases. This is due to an enhancement of the geomet- 
ric pinning of vortices and antivortices as R gets smaller. 
However, for 1.1 < R < 1.5 the behavior of ji(R) be- 
comes qualitatively different: it decreases with reducing 
R. This decrease is related to the fact that at those radii 
the distance between vortices, located at <p w 7r/2, and 
antivortices, located at <p w 37r/2, becomes comparable 
to the vortex size, so that their annihilation is facilitated 
by reducing R. At R < 1.1 even rather weak external 
current densities j e are sufficient for complete annihila- 
tion of vortex- antivortex pairs. In this case, the onset of 
a resistive state is determined by the nucleation of new 
pairs, which requires larger applied current densities j e 
at smaller R. As a result, in this range of R the critical 
current density rapidly increases with decreasing R. 

The calculated dependence of ji on the applied field 
Bo± [Fig. El(b)] is qualitatively similar to that found in 
Ref. [3 in the limit d — >• and for applied fields below H C 2 
(Bq± < 2). Increasing the shell thickness d leads to an 
increase in the critical current density ji due to partial 
screening of the applied field Bo± by supercurrents, so 
that a higher current must be applied in order to pro- 
duce a Lorentz force sufficient for vortex/antivortex de- 
pinning. The shape of the curves ji(Bo±) reveals several 
regimes, (i) At relatively low magnetic fields, no vortices 
appear in the cylinder, so that the values of j± actu- 
ally correspond to the transition from the Meissner state 
to the normal state, (ii) With increasing applied mag- 
netic field, the critical current density for nucleation of 
vortex- antivortex pairs, becomes smaller than the depair- 
ing current density so that j\ corresponds to the onset 
of vortex-antivortex nucleation, propagation, and anni- 
hilation, (iii) At even higher fields, vortices and antivor- 
tices are present in the cylinder already at j e < j\. The 
critical current density decreases monotonously with in- 
creasing Bq± [both due to an increase of the depinning 
Lorentz force, proportional to £?o_l, an d an enhanced mu- 
tual repulsion of the increasing number of vortices (an- 
tivortices)]. (iv) When approaching the second critical 
field (Bq± = 2), the number and density of vortices (an- 
tivortices) become so large that the distances between the 
vortex cores are smaller than the vortex size. The result- 
ing enhancement of mutual repulsion between vortices 
(antivortices) leads to a rather fast decrease of ji with 
Bo± in this field range. At Bo± = 2 the curves ji(Bo±) 
exhibit a pronounced minimum, (v) At Bq± > 2, the crit- 
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FIG. 2. (Color online) (a) Critical current density j± as a 
function of the cylinder radius R at L = 50, d —> 0, B ± =0.9, 
and Bon = 0. For comparison, the critical current density ji 
is shown also for a flat stripe with width 2ttR, length L = 50 
and thickness d — >> in a perpendicular magnetic field Bo± = 
0.9. (b) Critical current density j\ as a function of B ± at 
R — 3, L = 50, Bo\\ = 0, and different d. The dashed line 
roughly shows the boundary between the ranges of parameters 
where the critical current density ji corresponds to the onset 
of vortex-antivortex nucleation and recombination (below the 
line) or to the transition from the Meissner state to the normal 
state (above the line). 




3# n /2 20 



7' e =0.24 
■— B increases 
B Q decreases 



AAA/ 




FIG. 3. (Color online) Time- averaged voltage drop V between 
points 1 and 2 [see Fig. 1(c)] in a hollow cylinder with L — 75, 
R — 3, and d — >• as a function of an increasing magnetic field 
-B 1| for B ± = 0.35 and different applied current densities j e . 
Inset: Voltage drop V as a function of the applied parallel 
magnetic field in the cases of increasing (squares) and de- 
creasing (stars) field B \\ for L = 75, R = 3, d — » 0, and 
B ± = 0.35. 



ical current density first sharply increases and then slowly 
decreases. In this field range, normal regions are formed 
in the cylinder around <p = n/2 and <p = 37r/2, and the 
vortex dynamics involves the entrance of vortices and an- 
tivortices from the normal regions to the superconduct- 
ing regions, which requires relatively high applied cur- 
rent densities. Obviously, only the range of parameters 
corresponding to regimes (ii) to (iv) can be relevant for 
revealing the effect of persistent currents, induced by the 
field £> ||, on vortex dynamics. As illustrated by Fig. [2j 
this range is sufficiently wide. 

In Fig. [3] we plot the calculated time-averaged volt- 
age drop V between points 1 and 2 [see Fig. DJa)] as a 
function of an increasing parallel magnetic field B Q n for 
different values of the externally applied current density 
j e . The accuracy of the shown V- values is determined by 
a limited averaging time in the performed calculations. 
When switching B Q u to a new value, a time interval t r , 
ranging from 100 to 500, is reserved in the computational 
program for transient processes. Within this time inter- 
val, no calculation of V is performed. Then the voltage 
drop V is averaged over a time interval t ave r, ranging 
from 500 to 5000 in the present calculations. However, 



our analysis shows that in the multi- vortex system under 
consideration the full period of vortex-antivortex genera- 
tion/recombination processes can be rather long, so that 
in general this period is not much smaller than the used 
time intervals t aver . The error bars, shown in Fig. [3j cor- 
respond to the estimates, obtained using different values 

Of t r and taver- 

As seen from Fig. [3j there are well pronounced os- 
cillations of the calculated voltage V versus B \\. The 
oscillation period is typical for the AB oscillations and 
equals ^o/(ttR 2 ). The voltage V takes maximum val- 
ues for vanishing persistent currents induced by the field 
B ||, i.e. for ttR 2 B 011 w (2n + l)$ /2 (n e Z), while 
the sharp minima of V approximately correspond to the 
maximum magnitude of the persistent currents, i.e. to 
ttR 2 B \\ = (n + l/2)$ (n € Z). In other words, the 
shape of these oscillations is "inverted" as compared 
to the Little-Parks oscillations 13 . The oscillating be- 
havior of V(£? ||), shown in Fig. [3j can be explained 
in terms of the distortion of vortex/antivortex trajecto- 
ries due to the Lorentz forces caused by the persistent 
currents and the magnetic field Bq±. In Figs. SJa) to 




FIG. 4. (Color online) Panels (a) to (c): Distributions of the 
streaming parameter S for L = 75, R = 3, d — >• 0, B ± = 0.35, 
j e = 0.26, ti — t\ — 15 and different values of the applied 
parallel magnetic field: Bqu — (a), Bqu — 0.1 (b), and 
S 1| = 0.125 (c). Panel (d): Snapshot of the square modulus 
of the order parameter in the cylindrical shell with L = 75, 
R = 3, and d -> at S ± = 0.35, j e = 0.3, and B u = 0.2. 



IH[c), we plot the distributions of the streaming parame- 

r * 1 1/2 

ter S = (t 2 - h)- 1 J^(d\tlj\ 2 /dt) 2 dt , introduced in 

[2] to visualize the vortex/antivortex motion. As im- 
plied by Fig. SJa), in the case of B Q u = 0, and hence 
in the absence of the corresponding persistent currents, 
the vortex/antivortex trajectories in the middle part of 
the cylinder follow a circular cross-section of the cylin- 
der, with vortices moving from <fi = to tt and antivor- 
tices moving from (j) = 2ix to tt all at the same value 
of z. Vortex- vortex interactions cause only slight devi- 
ations from ideal circles. However, in the cases when 
the field B \\ induces strong diamagnetic [B \\ =0.1; see 
Fig. H(b)] or paramagnetic [B \\ = 0.125; see Fig. H(c)] 
currents in the cylindrical shell, these currents cause a 
rather pronounced deformation of vortex/antivortex tra- 
jectories, so that the path of a vortex (antivortex) from 
= (2tt) to (j) = 7r partly follows the ellipse traced out 
by a slanted cross-section of the cylinder. This represents 
an appreciably longer path than for B$\\ = 0. 

The behavior of V(Bq\\) at the applied current density 
j e = 0.24, which only slightly exceeds the critical current 
density j\ = 0.2392 at B$\\ = 0, implies that - in addi- 
tion to oscillations of V(B u) at a fixed j e -value - also 
the critical current density j\ is an oscillating function 
of the applied parallel magnetic field £> || =0. However, 
our calculations show that the corresponding oscillation 
amplitudes are very small. For the parameters under con- 
sideration the critical current at ttR 2 Bq\\ = (n + l/2)<l>o is 
ji = 0.2416, i.e. the oscillation amplitude for j\ does not 
exceed 0.6%. Thus we conclude that the voltage V(B Q u) 
is a better indicator than the critical current to reveal 
the AB oscillations in vortex trajectories. 



Besides pronounced oscillations of the voltage V(B \\) 
with period &o/(7rR 2 ), one can see in Fig. [3] several 
smaller features in the behavior of V(B Q n). Many of 
those features are irregular and reflect the limited compu- 
tation accuracy. However, some others are rather regular 
and periodic in B \\. In particular, for j\ = 0.24, small 
peaks of V(B Q u) emerge at magnetic fields slightly above 
B \\ = (n + l/2)^o/ (ttR 2 ). Those features can be at- 
tributed to the fact that the (time-averaged) number of 
vortex-antivortex pairs in the cylinder depends not only 
on the applied current density [see Figs. 0Jb) and [U(c)] 
but also on the parallel field £> ||. At the same time, the 
time-averaged voltage drop V as well as the critical cur- 
rent ji are, of course, sensitive to the number of vortices 
in the cylinder. This sensitivity leads, in particular, to 
the appearance of the peaks of V(-B ||) for je = 0.24 and 
S || slightly above (n + l/2)$ /(7r# 2 ). 

From Fig. O one can also observe that the curves 
V(B \\) are not fully symmetric with respect to the points 
5 || = n^/(uR 2 ) or B || = (71 + l/2)$ /(^R 2 ). Thus, 
in Fig. [3l the aforementioned weak peaks of V(B u) at 
ji =0.24 are present at magnetic fields just above B Q u = 
(n + 1/2)^>o/(tt^ 2 ), but no similar peaks appear below 
B \\ = (n -\- l/2)^o/ (ttR 2 ). In our numerical simulations, 
we ramp up the magnetic field B \\ and calculate how the 
order parameter adapts to the increased field. Geomet- 
ric pinning of vortices in the cylinder leads to hysteresis 
effects and results in asymmetries in the curves. A sim- 
ilar pinning in superconducting spherical nanoshells has 
been shown to cause a pronounced hysteresis in the de- 
pendence of the number of vortex-antivortex pairs on the 
applied magnetic field£. Also in the cylindrical shells un- 
der consideration the number of vortex-antivortex pairs 
and hence the curves V(B \\) demonstrate hysteretic be- 
havior. The presence of such a hysteretic behavior is 
illustrated in the inset of Fig. [3l where we compare the 
dependences V(-Bo||)? calculated for increasing and de- 
creasing fields B \\. As seen from this inset, in the case 
of a decreasing field £> ||, small peaks of V(B \\) appear 
to the left from B \\ = (3/2)$ /(7r# 2 ). Within the error 
bars, the whole pattern of V(Bq\\), where the field B \\ 
goes either up or down, looks symmetric with respect to 
Boy = (3/2)<J? /(7Ti? 2 ). 

As further seen from Fig. [3l the shape of the curve 
V(Bq\\) becomes more regular and symmetric when in- 
creasing the applied current density to the value j e = 
0.26, which is considerably higher than the critical cur- 
rent density j\. However, relatively pronounced ad- 
ditional features, caused by variations of the (time- 
averaged) number of vortex-antivortex pairs in the cylin- 
der, reappear in the curve V(B Q u) at j e = 0.3 (Fig. [3|), 
when - as illustrated by Fig. ffl^d) - some precursors of 
phase-slip-line formation can be already seen in the or- 
der parameter pattern. Remarkably, within the whole 
range of current densities considered here, the oscillation 
amplitude of V(Bq\\) remains appreciably large. 

From Fig. [3j the magnitude of the resistivity oscilla- 
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FIG. 5. (Color online) Time- averaged voltage drop V as a 
function of an increasing magnetic field £?o|| in the cases of 
(a) L = 30, R = 3, d = 0.272, B 0± = 0.6, j e = 0.28 and (b) 
L = 30, R = 6, d -> 0, B 0± = 0.2, j e = 0.2. 



follows from Fig. E(a), for L = 30, i? = 3, d = 0.272, 
ft = 0.77, Bq± = 0.6 and j e = 0.28 the oscillation 
magnitude is about 0.07. This magnitude tends to de- 
crease when increasing the radius of the cylinder and/or 
when decreasing its length down to values L < 2ttR [cp. 
Fig. Efb) to Fig. E]. Nevertheless, for L = 30, R = 6, 
d — »• 0, £?o± =0.2 and j e = 0.2 our calculations predict 
an oscillation magnitude as large as 0.01 [see Fig. Efb)]. 
Of course, in very short cylinders with R ~ 1, which 
cannot accommodate vortices, the predicted voltage os- 
cillations become impossible. 

It seems worth emphasizing that the material and geo- 
metric parameters, required to observe the predicted re- 
sistivity oscillations, are achievable experimentally. For 
example, with £ w 150 nm (close to the values of £(0) 
for Al hollow cylinders in Ref . 0) , the dimensionless pa- 
rameters R = 3 and d = 0.272 would correspond to a 
hollow cylinder with radius about 640 nm (a few times 
larger than the cylinder radii in Ref. [2|) and wall thickness 
about 60 nm (twice that in Ref. . 

IV. CONCLUSIONS 

To conclude, we have shown that in hollow supercon- 
ductor cylinders, subjected to a tilted magnetic field, 
the resistance, caused by vortex motion, should mani- 
fest measurable oscillations as a function of the magnetic 
field component parallel to the cylinder axis. This effect 
can provide a robust tool to probe experimentally the 
oscillating persistent currents, related to the Aharonov- 
Bohm effect, in a wide range of parameters, in particular, 
much below the superconducting critical temperature. 



tions is > 0.05 in the used units. This is about 0.5% 
of the normal-state resistivity, which equals 12 in our 
units^, so that the predicted oscillations should be ob- 
servable through 4-probe measurements similar to those 
reported in Ref. 0. Our calculations show that the self- 
inductance of the cylinders under consideration has a 
relatively weak effect on the resistivity oscillations; as 
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